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Abstract: The correlations between air temperatures, relative and absolute humidity, wind, cloudi- 
ness, precipitation and number of influenza cases have been extensively studied in the past. Be- 
cause, initially, COVID-19 cases were similar to influenza cases, researchers were prompted to look 
for similar relationships. The aim of the study is to identify the effects of changes in air temperature 
on the number of COVID-19 infections in Poland. The hypothesis under consideration concerns an 
increase in the number of COVID-19 cases as temperature decreases. The spatial heterogeneity of 
the relationship under study during the first year and a half of the COVID-19 pandemic in Polish 
counties is thus revealed. 


Keywords: COVID-19; cases; meteorological factors; temperatures; counties; Poland 


1. Introduction 


The impact of local meteorological factors on influenza morbidity has been observed 
in different geographical places in moderate northern latitudes. Correlations between air 
temperature, relative and absolute humidity, wind, cloudiness, and precipitation and the 
number of influenza cases have been extensively studied in various regions of the world 
[1-3]. Initially, the similarity between COVID-19 and influenza cases prompted research- 
ers to look for similar relationships. Mainly because some viral infections that occur dur- 
ing the cold seasons of the year are commonly called “colds”, which emphasizes their 
association with low temperature [4,5]. On the other hand, the risk of infection can para- 
doxically be reduced when the temperature drops to very low values [5-7]. Worldwide 
epidemiological analyses have shown a significant negative (weak to moderate) correla- 
tion (R) between temperature and the number of cases (incidence), mortality, recovered 
cases, and active cases of the COVID-19, with statistical significance (p = 0.0001-0.0030 for 
a 2-tailed test) and association strengths ranging from weak to moderate (R ranging from 
—0.205 to -0.332) [8]. 

Nevertheless, numerous other studies have found either positive or negative associ- 
ations of air temperature, humidity, and UV radiation with reported COVID-19 case num- 
bers [9-15]. Many of these studies were limited by short observation periods, either taking 
into account varying reporting biases, or failure related to the time lag between observed 
weather conditions and the acquired data reports concerning cases of infections and 
death. The reported observations were confirmed by laboratory evidence on the stability 
of SARS-CoV-2 as a function of temperature and humidity, proving that the virus’s half- 
life in human nasal mucus and sputum is shorter under higher temperature conditions 
and relative humidity than under conditions of lower temperature and relative humidity 
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[9]. Moreover, during winter people spend more time indoors, which may make the virus 
act unpredictably (either as isolation or growing factor, which reduces or raises the num- 
ber of infections). Some reports demonstrated that air temperature and relative humidity 
had a linear statistical relationship with the number of COVID-19 infections in different 
geographical regions in India [16]. On the other hand, a review of numerous published 
studies suggested that temperature is a driver of COVID-19 transmission intensification, 
but not necessarily the only one [17,18]. Some scientists anticipated that a warmer climate 
would not have any protective impact on current susceptibility to the virus, for example, 
in sub-Saharan Africa population centers [19]. There are also scientific papers, in which 
the results presented a correlation between the spread of droplets, and temperature be- 
cause respiratory viruses survive better from cold to dry weather [20]. The seasonality 
similarity of COVID-19 to other human coronaviruses and influenza let to assume hy- 
pothesis with a peak during winter months because of lower temperature, lower absolute 
humidity, and lower indoor relative humidity [18]. However, mixed evidence is available 
regarding the effect of changes in temperature on the number of cases of SARS-CoV-2 
infections [21]. Several papers reported an association between an increase in temperature 
and a decrease in newly infected COVID-19 cases [10,12,18,22-26]. The opposite relation 
has also been reported [26,27], as well as there being no relations at all [28-31]. Steiger [21] 
showed that weather, especially temperature (which has a reducing effect on case num- 
bers) and rainfall (which increases case numbers) affects the reported number of SARS- 
CoV-2 infections. 

This study aimed to investigate the spatial patterns of COVID-19 cases (SARS-CoV- 
2 infected, fatal, and recovered) and the effects of climate and bioclimate on observed 
changes in infected cases at the state and county levels in Poland. We hypothesized that 
with the increasing values of (1) air temperature, (2) dew point temperature, (3) ultraviolet 
radiation, (4) mean radiant temperature, and (5) Universal Thermal Climate Index, the 
COVID-19 growth rates decrease. This also prompts a hypothesis for pandemic seasonal- 
ity and the recurrence of pandemic waves. 


2. Materials and Methods 


The spatial units used in the study were the set of 380 Polish counties (powzat, the 
second-level unit of local government and administration, and NUTS-4 in the EU’s No- 
menclature of Territorial Units for Statistics). This is the most detailed level for regular 
daily reports on COVID-19 cases officially published in Poland. All the daily COVID-19 
cases (infected, fatal, and recovered) were aggregated into the weekly totals, retaining 
their spatial division. The observation period covered the 65 weeks from March 2020 till 
July 2021. There is a data gap only from the end of November to mid-December 2020 
(weeks 42-47) when counties lacked reported data. The main approach to spatial and time 
study is shown in Figure 1. 
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Figure 1. The study approach to the weekly relationships between COVID-19 cases and certain meteorological parameters 


(March 2020-July 2021). 
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The progress of the COVID-19 pandemic in Poland is described in Section 3. Two 
enormous waves of COVID-19 pandemic are revealed; the first: from October 2020 to De- 
cember 2020, and the second: from March 2021 to April 2021. Data at the national level 
were additionally acquired using the CovsirPhy Python library and a dedicated applica- 
tion (Python script). The application used downloaded data from the COVID-19 Data Hub 
[32-34]. A cross-check of the data was validated using officially published data from the 
Polish Ministry of Health [35]. All daily reported data concerning the number of COVID- 
19 cases —SARS-CoV-2 infections cases, failed cases (deaths), and recovered — were aggre- 
gated into weeks. 

The meteorological and biometeorological parameters were obtained from the 
ECMWE ERA5 and ERA5-HEAT reanalyzes [36,37]. The reanalyzes are hosted at the Cli- 
mate Data Store (CDS: https://cds.climate.copernicus.eu/, access: 20 July 2021), imple- 
mented by the European Centre for Medium-Range Weather Forecasts (ECMWF) as part 
of the Copernicus Climate Change Service (C35). The data are currently available from 
1979 (but will be pushed back to 1950) to the present with a frequency of up to 1 handa 
spatial resolution of 0.25° x 0.25° (approximately 28 km x 28 km). The data are stored in 
GRIB (GRIdded Binary) format, a World Meteorological Organization (WMO) standard 
for the storage and distribution of gridded data. 

The variables used in the study were 2 m air temperature (Ta), 2 m dew point tem- 
perature (Td), downward UV radiation at the surface (UV), mean radiant temperature 
(MRT), and Universal Thermal Climate Index (UTCI) were used. The MRT is defined as 
the “uniform temperature of a fictive black-body radiation enclosure (emission coefficient 
€ = 1) which would result in the same net radiation energy exchange with the subject as 
the actual, more complex radiation environment” [38] for an individual placed in a given 
environment in a given position and clothing. It is converted into a temperature using the 


Stefan—Boltzmann Law. 
T ca [BS 273.15 1 
mrt Epô z ( ) 


where: Tint [°C], Ssr—absorbed flux density, ô—the Stefan-Boltzmann constant, and €p- 
emissivity of the clothed human body [38]. MRT is widely applied in various disciplines, 
from engineering, urban studies, and climate change to physiology and public health [39]. 
In biometeorology, it is an input variable for thermal comfort and stress indices, parallel 
to UTCI, which is also used in this study. The UTCI is defined as the reference condition’s 
air temperature (Ta) that causes the same model response as actual conditions. The devi- 
ation of UTCI from the air temperature depends on the actual values of the air tempera- 
ture, MRT, wind speed and humidity [40]. Weekly means of all variables were calculated 
using the ISO 8601 calendar (from the tenth week of 2020 to the thirtieth week of 2021). 

The steps taken in preparing the data involved the recalculation of input data (di- 
vided into weeks) to convert all the temperature data layers for Poland, expressed in Kel- 
vin, into Celsius degrees. Next, using the zonal statistics function in GIS software, all ras- 
ter weekly temperature data (in degrees Celsius) were aggregated into the attribute tables 
using a template for a vector map of Polish counties, cutting off the data outside the state’s 
boundaries. The general approach to the data analysis involved comparing of the number 
of infected people, aggregated by weeks and by county in Poland (from week 10 of 2020 
to week 30 of 2021), with the maximum, mean, and minimum temperatures and volumes 
of radiation observed (see section 4). Additionally, the studied relationships were evalu- 
ated spatially over the observed time period, calculating values for the Pearson’s correla- 
tion coefficient separately for each county over time (see section 5). 

The main approach to the spatial and temporal study is presented in Figure 1. In the 
first stage, the general approach involved the general analysis of the summarized and 
mean values figures for the whole of Poland, and, later, grouped into weeks and counties, 
aimed to identify the most correlated independent meteorological variable suspected of 
having an impact on the rise in the number of infections. The second, detailed stage of the 


Remote Sens. 2021, 13, 4946 


4 of 16 


analysis a concerned the spatial and temporal analysis of the selected relationships. Pear- 
son’s correlation coefficients (R) and significance, checked using 2-tailed p-values over 
time, were computed using SPSS (Statistical Package for the Social Sciences). The GeoDa 
package was used for spatial analyses. 

The steps in the study involved the following: (i) an analysis of the general statistics 
for the whole of Poland concerning COVID-19 and meteorological factors (aimed at iden- 
tifying any significant effects of meteorological conditions on COVID-19 cases); and (ii) a 
detailed spatial and temporal statistical analysis of the selected meteorological factors, in- 
cluding the map algebra zonal statistics tools, which enabled these data to be linked with 
COVID-19 infections in the form of data cube attributes grouped into weeks and counties, 
and cartographical presentations over time. Map algebra is an appellation for a collection 
of tools introduced at the end of the twentieth century [41] for raster processing on several 
collections of images. Map algebra tools processing can be used for physical [42] and so- 
cial-economical [43] studies for various spatial sizes in order to create local, focal, neigh- 
borhood, zonal, or global statistics [44]. 

Next, the spatial diversity of the relationship between the number of infections and 
the mean radiant temperatures were examined as a Bivariate Local Moran I (Index) by 
selecting particular, certain weeks when the number of infections rose globally in Poland, 
i.e., the peaks of COVID-19 pandemic waves. These analyses aimed to distinguish four 
types of regions (counties) in which spatial lags between mean radiant temperatures and 
infections occurred in certain weeks (High-High, Low-Low, High-Low, Low-High) and 
their spatial heterogeneity. 

The set of computer tools and computer applications used in the study involved GIS 
(ArcGIS, QGIS, and GeoDa), Python scripts (aimed transformation of meteorological 
data), and statistical software (SPSS). 


3. The COVID-19 Pandemic in Poland 


The first infection was reported on 4 March 2020. Two early moderate waves (at the 
beginning of 2020 when compared to later ones) were reported, followed by two main 
significant waves of infections, one in autumn 2020 and the other in spring 2021 (Figures 
1 and 2); this despite the efforts of national healthcare services, including regional and 
national lockdowns, organizing both a stationary and field hospital network, and starting 
the vaccination of the population from January 2021. The progress of the COVID-19 pan- 
demic in Poland is presented in Figures 2 and 3. 
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Figure 2. Poland, COVID-19 cases over time with monotonic increasing complemented recovered data (as of June 04, 2021, 
results: CovsirPhy Python Application). 
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Figure 3. Stacked, cumulative number of infections by county in Poland per week (from March 2020 to July 2021); data 
source: Ministry of Health, Poland (counties identified by statistical numbers, four-digit TERYT string [35]). 


4. Climatic and Bioclimatic Conditions in Poland 


According to the Köppen-Geiger climate classification, almost the whole territory of 
Poland is classified as having a humid continental climate with warm summers (Dfb) [45]. 
Only a small area in the high mountains (on the southern and southwestern borders) is 
classified as having a subarctic climate (Dfc). 
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4.1. Temperatures 


According to the 2020 annual report by the Institute of Meteorology and Water Man- 
agement— National Research Institute thermal conditions for all months (except May and 
July 2020) may be classified as at least slightly warm across all physiogeographic regions 
of Poland. February and August were extremely warm months almost everywhere, May 
was extremely cool, and July was thermally normal in comparison to the empirical distri- 
bution of average temperatures for particular months, seasons, and years since from 1981 
to 2010 (Table 1). 


Table 1. Global average monthly air temperatures anomalies (1981-2010), and classification of thermal conditions in Po- 
land by month in 2020 (source: Institute of Meteorology and Water Management, National Research Institute, 2021) [46]. 


JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC YEAR 
2020 22 40 45 87 110 17.7 182 198 151 105 56 19 99 
Anomalies, 1981-2010 37 46 17 06 -2.3 1.7 -01 20 18 19 22 22 1.6 


Detailed monthly temperatures over the observed time were obtained by estimating 
mean temperatures based on values of (weekly) raster layers acquired from the Climate 
Date Store (NetCDF format), estimating the spatial extent of Poland’s state boundaries 
with wider margins (west longitude = 13°00’E; north latitude = 55%54'5”N; east longitude 
= 25°00’E south latitude = 48°05'5” N). 

The meteorological (temperature) variables taken into account in the study hereafter 
are the near surface air temperature (nsat) and the dew point temperature (dpt). The bio- 
meteorological parameter taken as the variable is the mean radiant temperature (MRT). 

The range of the weekly minimum, mean, and maximum values for 2 m air temper- 
ature (nsat), 2 m dew point temperature (dpt), and mean radiant temperature (rt) over the 
observed time, globally in Poland, are presented in Figure 4. The weekly minimum, mean 
and maximum values reflected typical seasonality for temperate latitudes. 


Minimum, mean and maximum temperatures. by week in Poland 


(January 2020—July 2021, t: radiant temperature, nsat: near surface air temperature, dpt: 2m dew point temperature) 
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Figure 4. Weekly minimum, mean, and maximum temperatures in Poland (January 2020—July 2021) [°C]; descriptions: rt- 
radiant temperatures (MRT); nsat-2 m near surface air temperatures; dpt-2 m dew point temperatures. 
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These data have been recurrently and iteratively grouped by county using the (GIS) 
zonal statistics tool with a zone (template) map of Polish counties. The results were a set 
of mean temperatures by county and by week over the observed period. These data were 
confronted with the number of SARS-CoV-2 infected cases grouped by county and by 
week. 


4.2. UV Radiation at the Surface 


The range of the weekly minimum, mean, and maximum values for the volume of 
ultraviolet radiation (J per sq. m) over the observed time, globally in Poland, are presented 
in Figure 5. UV radiation was taken into account in the study due to its above-indicated 
importance in the literature and its potential impact on diminishing the spread of corona- 
virus SARS-CoV-2. UV minimum, mean, and maximum values gradually increase from 
winter to summer with the inverse occurring from summer to winter. 


4.3. Universal Thermal Climate Index 


The only composite measure—UTCI—was also examined because it is one of the pri- 
mary bioclimatological indicators. Due to the possibility of not being able to determine if 
a single meteorological factor influences the spread of coronavirus SARS-CoV-2, this com- 
plex measure also should be studied. The range of the weekly, minimum, mean, and max- 
imum values for UTCI over the observed time, globally in Poland, are presented in Figure 
6. 


Downward UV radiation at the surface in Poland by weeks 


January 2020-July 2021 


— mean LA — minty — maU 





Figure 5. Downward ultraviolet radiation at the surface (UV) [J/m2]: weekly minimum (min), mean, and maximum (max) 
values in Poland (January 2020-July 2021). 
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Universal Thermal Climate Index UTCI by weeks in Poland 
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Figure 6. Universal Thermal Climate Index (UTCI) [°C]: weekly minimum (min), mean, and maximum (max) values in 
Poland (January 2020-July 2021). 


5. General Approach 


The general approach to the study involved estimating the Pearson correlation coef- 
ficients for pairs of variables: between each of the chosen meteorological and biometeor- 
ological factors (variables) described above, and the number of SARS-CoV-2 infections by 
week and by county over the examined time. These procedures were aimed at identifying 
relationships and showing the global relationship within Poland. The general results con- 
firmed the primary hypothesis that the number of infections decreases as the tempera- 
tures, UTCI, and UV radiation grow (Table 2, Figure 7A,B). 


Table 2. Pearson’s correlation coefficients for temperatures, UV radiation, UTCI, and total number of SARS-CoV-2 infec- 
tions and failed cases (deaths) in Poland over time in weeks (March 2020—July 2021). 


Correlations 
In- Fail mea min max mean min max mean min max mean min max mean min max 
fected ed nrt rt rt nsat nsat nsat dpt dpt dpt UTCI UTCI UTCI UV UV UV 


Pearson’s 1 0.841 -0.694 -0.703 -0.676 -0.655 -0.662 -0.614 -0.491 -0.527 -0.521 0.651 ** -0.636 -0.660 -0.666 -0.651 -0.657 
Colao PER m m m e se se se A se se se 
In- - 
fected ner 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
N (weeks) 65 65 65 65 65 64 64 64 64 64 64 65 65 65 64 64 64 
Pearson 0841 1 —0.667 -0.688-0.630 -0.644 -0.643 -0.585 -0.498 -0.520 -0.510 0.649 ** -0.625 -0.636 -0.585 -0.575 -0.586 
Cancion ER se m se e A se se A e se A 
Failed Sig. (2- 
tailed) 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 


N (weeks) 65 65 65 65 65 64 64 64 64 64 64 6 65 65 64 64 64 
** Correlation is significant at the 0.01 level (2-tailed); Descriptions: MRT rt—radiant temperatures, nsat-near surface 
air temperature, dpt-2 m dew point temperatures, UTCI-Universal Thermal Climate Index, UV-Downward UV radi- 
ation at the surface. 
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Figure 7. (A). The general relationships between mean radiant temperatures and number of SARS-CoV-2 infections, ag- 
gregated in weeks in Poland —3D scatterplot (March 2020-July 2021); (B). Scatterplot matrix: the relationships between 


mean radiant temperature (MRT) and number of SARS-CoV-2 infections aggregated in weeks over time. 


The complex chart (Figure 8) presents the number of COVID-19 cases related to av- 
erage week mean radiant temperatures (MRT) for each county, for each week during the 
observed period. There is a data gap only from the end of November to mid-December 
2020 (weeks 42-47) when the counties lacked reported data. The chart reveals evident sea- 
sonality with peaks in winter and early spring (not including March 2020, the starting 
time of pandemic diffusion). 
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Figure 8. Poland. COVID-19 infections for weeks and by counties vs weekly mean radiant temperatures (from March 


2020, week 11 to July 2021, week 30). 
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The spatial heterogeneity of the relationships between COVID-19 cases and the 
weekly mean radiant temperature during the whole observed time is revealed on the map 
presenting the set of estimated Pearson's coefficients for correlations grouped by county 
(Figure 9A). The separate correlation coefficients (R), as well as determination coefficients 
(R2), were calculated for each county over the observed time. They are also presented on 
the map (Figure 9B). 
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Figure 9. Poland. (A). The Pearson’s correlation coefficients for weekly mean radiant temperature and total number of 
COVID-19 infections grouped by county over time. (B). (R2) coefficients of determination for weekly mean radiant tem- 
perature and total number of COVID-19 infections grouped by county over time (March 2020-July 2021). 


6. Results 


It is important to note that the observed period begins at week 10 of 2020 and ends at 
week 30 of 2021 (March 2020 to June 2021). 


Results of Detailed Spatial Analyses 


Correlation coefficients between the weekly mean radiant temperature and number 
of COVID-19 weekly cases by county, reveal some northern, eastern, and western regions 
are where MRT seemed to matter (Figure 9). Analysis of the map of coefficients of deter- 
mination (squares of correlation coefficients, Figure 9B) by county, and analysis of spatial 
autocorrelation, confirms the observed relation over time, grouped by county (Univariate 
Global Moran I = 0.588). However, there are no lower or upper outliers. 

Examination of spatial autocorrelation and statistical significance, using local bivari- 
ate Moran’s I (BiLISA) coefficients, assumed the mean radiant temperature to be the inde- 
pendent variable, and the number of infections (Summed by weeks) to be the dependent 
variable (Figure 10). Finally, K-means clustering shows the spatial grouping of counties 
assuming categorization into five unique groups of counties based on the relationship of 
both observed variables: weekly mean radiant temperature and the sum of weekly num- 
ber of infected (COVID-19 cases, Figure 10). 
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Figure 10. Relationship between weekly mean radiation temperature and COVID-19 cases (SARS-CoV-2 infections) by 
county, over selected weeks in Poland. (A). Bivariate Local Moran’s I (BiLISA) cluster maps by selected weeks. (B). Biva- 
riate Local Moran’s I (BiLISA) significance maps by county over selected weeks. (C). Additional statistical information. 
(D). K-means clusters (mean radiant temperature, number of COVID-19 cases) over selected weeks. 


7. Discussion and Conclusions 


The nexus of factors influencing the spread of the SARS-CoV-2 pandemic is so com- 
plex that the identification of (some) determining factors of COVID-19 spatial diffusion is 
significantly hampered. It has been shown that COVID-19 characterizes of specific dy- 
namics and is far more transmissible [47]. The mobility of people and their number, as 
well as the density of the population and, after all, the intensity of direct social contact, in 
turn, results in an increase in spatial interactions and the number of infections [37]. Spatial 
natural conditions (inter alia temperatures) are also recognized as primary factors in en- 
vironmental risk. Early in 2020, researchers found that all cities experiencing significant 
outbreaks of COVID-19 had very similar winter climates, with an average temperatures 
of 5 to 11 degrees Celsius, an average humidity levels of 47 to 79 percent, and a narrow 
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east-west distribution along the same 30-50 N latitude [48]. Most of studies suggested that 
climate may only modulate the transmission of severe acute respiratory syndrome coro- 
navirus 2 (SARS-CoV-2) and it remains unclear whether seasonal and geographic varia- 
tions in climate can substantially alter the pandemic’s trajectory [49]. 

The general approach to the statistical analysis of the relationship between mean ra- 
diant temperature and number of SARS-Cov-2 infections over time in Poland allows us to 
infer the need for further research regarding the hypothesis of seasonality of the COVID- 
19 pandemic. 

The research showed a strong correlation between mean radiant temperature and the 
number of infections over time in Poland. All the examined meteorological and biocli- 
matic factors, especially temperature, were inversely correlated to the total number of in- 
fections in Poland over time (for a sequence of 65 weeks, from March 2020 to June 2021). 
The Pearson’s correlation coefficient value (R = —0.7) indicates a decrease in infections as 
mean radiant temperatures rise over time, which is the meteorological factor most 
strongly, inversely correlated to infections. Detailed statistical analysis for particular 
counties, i.e., the estimation the analogous relationships over time (dependent variable— 
infections over selected examined weeks), allows us to conclude that the northern regions 
of Poland, as well as some eastern counties, showed similar strong inverse correlation 
values, but that they are more scarcely populated than central and southern Poland. How- 
ever, one should emphasize that the obtained results require interpretation within the 
context of the climatic, demographic, and topographic conditions characteristic for the 
territory of Poland and the analyzed period. Making inferences about the periodicity of 
the COVID-19 pandemic requires continued and more extensive research, both in terms 
of time (long-term) and area (a continent-wide region). 

The observed relationship between MRT and infections revealed in the general ap- 
proach, and for the northern and eastern counties of Poland over time, was unconfirmed 
during the study of bivariate detailed spatial diversity by county for particular weeks. The 
spatial diversity of MRT in Poland is uncorrelated with infections within the set of coun- 
ties for any selected week. There had to be other factors controlling the outbreaks of pan- 
demic during short periods of time (week) [47]. However, observations of the spatial ar- 
rangement of the first two clusters at the beginning of pandemic and the last weeks of 
examined time (reversal of pandemic) are interesting, and showed a strict regionalization 
and spatial consistency. This can be explained numerically: the number of infected people 
was relatively small in the prevalent number of counties, but only a few pandemic out- 
breaks were revealed and mostly in or around urban areas and cities, as well near state 
border crossings. 

It has to be remembered that the vaccination program against COVID-19 started in 
Poland on 27 December 2020 (week 52). However, the effectiveness of these vaccines at 
preventing infection with new virus variants, which have emerged in the last year, is not 
fully defined. Moreover, data are scarce on the risk of Delta variant transmission by vac- 
cinated individuals with mild infections. This shows that the factors influencing the 
spread of the COVID-19 pandemic may be more complex. Also, there is information that 
vaccination does not protect completely against infection, but is mostly effective at pre- 
venting severe disease and deaths from COVID-19. Vaccination is the best way to control 
the pandemic. There are now 99 vaccines in clinical trials, 25 have reached phase III effi- 
cacy studies and 18 have received some form of approval for use (as of July 2021) [50]. 
Despite authorization having been granted for multiple vaccines, as the ongoing global 
outbreaks demonstrate, the pandemic is far from over. SARS-CoV-2 continues to evolve 
under immune selective pressure, and while transmission levels remain high, there is the 
increased likelihood of vaccine escape variants evolving [50]. Coronavirus SARS-CoV-2 
has proven to be much more prone to mutations than researchers first thought, and more 
variants are emerging all the time. Six months is not much time to collect data on how 
durable vaccine responses will be. Recent research has shown that susceptibility to infec- 
tion with Delta, a highly transmissible variant, increases with time, as early as 2-3 months 
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after vaccination, which is consistent with observed waning protective immunity [51]. 
This indicates that the earliest groups to be vaccinated are more likely to have less robust 
immune responses to infection. Key clinical trials for currently authorized vaccines have 
determined whether the inoculations could safely avert symptomatic disease in individ- 
uals [52]. Studies have suggested that although vaccination reduces the risk of Delta var- 
iant infection and accelerates viral clearance, fully vaccinated individuals who become 
SARS-CoV-2-positive have peak viral loads similar to unvaccinated cases and can effi- 
ciently transmit infection [51]. Recent studies have also shown that there are no significant 
differences in cycle threshold values between vaccinated and unvaccinated, or asympto- 
matic and symptomatic groups infected with the SARS-CoV-2 Delta variant [53,54]. 

The epidemiological trends observed worldwide suggest a seasonal effect for the dis- 
ease that is not directly controlled by the genomic makeup of the virus [8]. Some research- 
ers suggest that air drying capacity (ADC) and ultraviolet radiation (UV) are probable 
environmental determinants in shaping the transmission of COVID-19 at a seasonal time 
scale [55]. Recent global studies have revealed consistent negative effects for both temper- 
ature and absolute humidity at large spatial scales for COVID-19 cases, but at finer spatial 
resolutions, these connections have been substantiated during the seasonal rise and fall of 
COVID-19 [56]. COVID-19 seasonality is more pronounced at higher latitudes, where 
larger seasonal amplitudes of environmental indicators are observed [57]. The general ap- 
proach of our study strongly suggests this hypothesis (see Table 2) in some, more prone 
regions in northern and eastern regions of Poland over the observed time (Figure 9). De- 
tailed analyses by week and by county (Figure 10) showed insignificant spatial relations 
over shorter intervals. The main social and demographic drivers of outbreaks in Poland 
during the first year of the COVID-19 pandemic were confirmed and described in a sepa- 
rate paper [47]. Secondly, “Weather may marginally affect COVID-19 dynamics (...) Fu- 
ture scientific work on this politically-fraught topic needs a more careful approach. ” [58]. 

This research has shown that one of the independent variables influencing the re- 
turning seasonal waves of the pandemic in Poland is the bioclimatic condition variable, 
particularly temperature (MRT). One should stress that the obtained results (Figure 9) in- 
dicate a strong correlation between biometeorological and meteorological parameters, 
and the number of COVID-19 cases in Poland. Simultaneously, they point to the spatial 
differentiation of this phenomenon: this correlation is much stronger in the north of Po- 
land, where the coefficient of determination reaches values as high as 80%. Thus, climate- 
related independent variables can account for a significant portion of the dependent var- 
iable. A time series analysis for the southern part of the country, characterized by a higher 
level of occupational mobility of inhabitants, should also consider, for example, data on 
mobility. Developing such a holistic multi-variable regression model will be the subject of 
further research. 
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